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SUMMARY 


^ An equation Is presented for calculating the heat flow required 
from the surface of an Internally heated windshield In order to 
prevent the formation of ice accretions during flight in specified 
icing conditions. To ascertain the validity of the equation,, 
conqsarlson is made between calculated values of the heat required 
and measured values obtained for test windshields in actxial flights 
in icing conditions. 

The test windshields were internally heated and provided data 
applicable to two common types of windshield configurations j namely 
the V-type and the type installed flush with the fuselage contours. 
These windshields were Installed on a twin-engine cargo airplane anil 
the icing flights were conducted over a large area of the 
United States during the winters of 1945-46 and 1946-47, In 
addition to the internally heated windshield investigation, some 
test data were obtained for a windshield ice -prevent! on system in 
which heated air was discharged into the windshield boundaiy layer. 

The general conclusions resulting from this investigation are 
as follows; 

1, The eunoxmt of heat required for the prevention of ice 
accretions on both flush- and V-type windshields during flight in 
specified icing conditions can be calculated with a degree of 
accuracy suitable for design purposes. 

2. A heat flow of 2000 to 2500 Btu per hour per square foot 
is inquired for conplete and continuous protection of a V-type 
windshield in fli^t at speeds up to 300 miles per hour in a 


2 


NACA TW Ho. 1434 


moderate cuBBtlotis icing condition. .For the. same degree of . protection 
and the same speed range, a value of 1000 Btu per hour per square 
foot suffices in a moderate stratus icing. condition. 

3 . A heat supply of 1000 Btu per hour per square foot is 
adequate for a flush windshield located well aft of the fuselage 
stagnation region, at speeds up to 3 OO miles per hoxir, for flight 
in both stratus and moderate cumuloue icing conditions. 

4. ' The external air discharge system of windshield thermal 
ice prevention is thermally inefficient and requires a. heat supply 
approximately 20 times that required for an internal system having 
the same performance. 


INTRODUCTION 

For several years the NACA has engaged in a broad research 
program on the problem of the prevention of ice formations on 
airplanes. Particular attention has been given to the utilization 
(in various thermal ice— prevention systems) of the available waste 
heat in the engine ejdiaust gases. . • 

One part of this icing research program has been concerned with 
the investigation of thermal means of windshield ice— prevention. 

The first satisfactory solution developed was the double— panel— type 
system described in reference 1 and tested on the Lockheed 12— A, 
Consolidated B-24, and Curtiss-Nright C-46 airplanes (references 1, 

2, and 3) • 

The tests of' reference 1 resulted in the tentative specifica- 
tion of a hea't— flow requirement of 1000 Btu per square foot per hour 
through the windshield outer surface. This value was based on data 
obtained for a T-type windshield at flight speeds up to I 50 miles 
per hour. The flight investigation of reference 2 was also conducted 
with a windshield of this configuration, and the specification 
appeared adequate. The heated— air flush windshield installation in 
the C— 46 airplane (reference 3) provided satisfactory protection, 
but did not serve as a check on the validity of the 1000 Btu per 
hour per square foot specification for all types of windshields^ 
because the windshield configuration had appreciably different icing 
characteristics than the V— types previously tested. 

Althoiigh these initial researches provided some information on 
windshield heat requirements, the results were empirical in nature 
and could not serve as a fimdamental basis for the prediction of the 
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heating requirements for windshield configurations and flight 
conditions different than those investigated. Accordingly, a 
fundamental windshield icing research was undertaken. This research 
consisted of an analytical study of the heating requirements for 
internally heated windshields, and subsequent flight tests to 
determine the applicability of the design "^equations resulting from 
the analysis. In the flight tests, the heat flow from the surfaces , 
of various special test windshields of different configurations 
could be closely controlled and measured, .Measurements were also 
made of the windshield surface temperatvires, boundary— layer 
profiles, flight conditions, and meteorological factors dviring 
icing in order to provide data on all the basic factors which could 
be expected to affect the heating requirement. 

In addition to the research directed toward the establishment 
of the heat— flow requirement in the case where the heating was 
supplied internally, a- secondary investigation conducted at the 
same time was concerned with the practicability of ice prevention 
by the means of discharging a jet of heated air into the windshield 
boundary layer. This device was initially installed in the 
C— 46 airplane (reference 3 ) as a means of augmenting the internal, 
double-panel windshield system with available primary air from the 
exhaust-gas heat exchangers. The test results with this Initial 
installation were sufficiently promising to warrant ftarther : 
investigation. 

An analytical approach to the external— discharge type of 
windshield thermal ice— prevention system was attempted, but the 
unknowns involved, 'such as the mixing of the 'heated Jet with the 
. boundary— layer air, precluded a reasonable prediction of the action 
of this system. Actual tests of external-discharge systems in 
icing conditions were considered necessary. It was hoped that 
sufficient data concerning the influence- of the veirious pertinent 
factors, such as jet flow rate, size, and temperatiire in known 
icing conditions, could be obtained to provide an indication of 
the practicability of the system and .possibly form the basis for 
empirical design equations. Controlled heated-air external-discharge 
systems, therefore, were investigated for both a V— type windshield, 
and a windshield which was flush with the -fuselage contours. 

The flight tests were conducted in clear 'air and in. natural 
icing conditions with a C— 46 airplane. In addition to the windshield 
research the flight tests included the detennlnation of wing and 
propeller heating requirements and a study of the meteorological 
factors conducive to icing. To obtain test data in natural icing 
conditions, for all parts of the research pro^am the C— 46 airplane 
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was operated "by the Ames Aeronautical Laboratory, Moffett Field, 
Calif., during the winters 1945-46 and 1946-47. 

The appreciation of the WACA is extended to United Air Lines, 
Inc., the United States Weather Bureau, and to the Air Ifeterial 
Command of the Army Air Forces for aid and cooperation in the 
research. In particular, the services of Jfajor James Murray, 

Army Air Forces, and Captain Carl M. Christenson and First Officer 
Lyle W. Eeynolds, United Air Lines, who served as pilots of the 
research airplane, were a valuable aid ,to the conduct of the 
investigation. • 


SYMBOLS ; 

The following nomenclature is vised throvighout this report: 
radius of water droplet, cenijimeters 
surface area, square feet • 
projected area, square feet 

projected height of windshield or half of a width of 
ribbon, centimeters 

specific heat of air, Btu per pound, degree Fahrenheit 
specific heat of water, Btu per povind, degree Fahrenheit 
radius of a sphere, feet 

a significant dimension of windshield used in determining 
Eeynolds number, feet 

e water vapor pressure, inches of mercury 

g acceleration due to gravity, feet per second, second 

h convective heat— transfer coefficient throvigh the wind- 

shield bovindary layer, Btu per hour, square foot, degree 
Fahreidieit 

J mechanical equivalent of heat, foot-pounds per Btu 

thermal conductivity, Btu per hour, square foot, 
degree Fahrenheit per foot 
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K 


I 

Ls 


m 

M 


P 

Pr 


<1 

Q 

r 

E 

s 

t 





T 

u 

U 

w 


dimensionless quantity oTstained frcan reference It- and 
defined in the symbols of this report as 

length of a windshield panel or a flat plate in the 
direction of local aif flow, feet 

latent heat of evaporation at surface temperature, Btu 
per pound of water 

liquid water content of the air, grams per cubic meter 

weight rate of water impingement per unit area, pounds 
per hour, square foot 

beirometric pressure, inches of mercury 

, nondimens ional 


I^andtls number 


k 


V 





unit rate of heat flow, Btu per hour> square foot 

rate of heat supply, Btu per hour 

X ys 

recovery factor equal to Pr ' for turbulent flow 
/ UrD \ 

Eeynolds number ( nondimens ional 

^ I - 

distance from the region of air stagnation, feet 
temperature, degrees Fahrenheit 

kinetic temperature rise of air, degrees Fahrenheit 


kinetic temperatvire rise of the water droplets, degrees 
Fahrenheit 

temperature, degrees Fahrenheit absolute 

local velocity inside boundary layer, feet per; second 

velocity, feet per 'second 

weight rate of water impingement, pounds per hour 
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W 

X 


X 

y 

a 
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s 

n 

9 


weight of heated air supplied, pounds per hour 

distance along windshield in the direction of air flow, 
inches 


0. 622 Ls / es-eo \ 

evaporation factor, defined as 1 + r ^ J 

distance normal to the windshield surface, inches 

angle between the plane of a windshield and a plane 
normal to the streamlines around the forebody of the 
windshield, degrees ' ; 

specific weight of air, pounds per cubic foot or grams 
per cubic centimenter 


specific weight of the water droplets, grams per cubic 
centimenter • ' • ■ 

emissivity of windshield panel 

efficiency of water impingement, percent 

hq lf the central angle of the total area of impingement 
on a spherical surface, degrees 


viscosity of air, gram seconds per square centimeter 
dimensionless quantity obtained from reference 4 and 


defined in the symbols of this report as 


iSy^Upb \ 

J 


Subscripts 

c calculated 

E experimental 

X local conditions Just outside the boundary layer 

o reference to ambient or free-stream air conditions 

s reference to windshield external surface conditions 

av 


average conditions 
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ANALYSIS 

Ice formations on aircraft occtir duing flight through formations 
of- supercooled droplets of water. Upon contact with an unheated 
surface, the water freezes and ice accumulates at a rate depending 
upon the efficiency with which the s\irface intercepts water droplets 
from the free air stream. With a heated surface the majority of 
the water evaporated is evaporated at the surface, although a few 
small drops may he evaporated in the boundary layer. The max imum 
rate of evaporation occurs when the heated surface is fully wetted, 
and it is for this condition that the following analysis applies. 

The unit heat flow from the outer surface of a windshield 
during flight in icing conditions can he considered as the sum of. 

losses, or 

(i = ^ + <ls ^ ( 1 ) 


due to forced convection 
due to Bvapoiation of the impinging water 
due to warming of the inqpinging water 
due to radiation to the sxirrounding atmosphere 
Each of these individual heat flows will he analyzed. 


four individual 


where 




heat 

loss 


heat 

loss 


heat 

loss 

<34 

heat 

loss 


Heat Loss Due to Convection ■ 

The eqmtion for the heat loss due to convection, including 
the effect of kinetic heating, is written 

q = h(t3-to-^tj, ) (2) 

£l 


where 




•_ Uq^-Ui^ 
. 2gJc^ 


.+ 




( 3 ) 
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Equation (3) is obtained from references 5 or 6 which also state 
that the recovery factor r equals Er^/® for turbulent, flcn^r. 

The first term of equation ( 3 ) represents the adiabatic 
heating or cooling of the air Just outside the boTindary layer 
caused by the change in velocity from Uq to Ui, The second 

term of the equation represents the heating of the air which occiurs, 
as the velocity is reduced from Ui to zero. Equation (3) can be 
approximated by the simplor form 

At^^ = 0.832 rfHo-Y ( 4 ) 

% VlOO/ 

Several equations for calculating the coefficient of convective 
heat transfer h from a flat' plate have been developed, and a 
comparison of these equations has been made in reference ?• For a 
turbulent boiindary layer, the general form of the coefficient was 
presented as ' • 

Ux7c.„ 

h = 0.0296 X 3600 ( 5 ) 

Following . the presentation of equation ( 5 ) in reference 7 j various 
specific forms of the equation are derived by expressing the 
properties of air as a function of the average tempearekture of the 
air in the boundary layer. The coefficient at any distance s from 
the leading edge of a flat plate, for the region of turbulent flow, 
is given as , 

h.0.51 ( 6 ) 

where 

*^av average temperature in the boundary layer 

Ui velocity outside the boundary layer at distance s from 
the. leading edge 

For a turbulent boundary layer extending from s=6 to 
a=l on a flat plate, the average coefficient is presented as 
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Heat Loss Due to Evaporation 

As mentioned previously, the maximum rate of evaporation 
occurs when the surface is fiilly wetted. Since this condition 
requires the maximum heat supply to maintain a given surface 
temperatvire, it is chosen for design. In order to establish the 
minimum heat requirement to prevent freezing of the water on the 
windshield, it is assumed that svtfficient heat will be supplied to 
maintain the temperatiure of the windshield sm^face and the water 
thereon at 32° F. Since water runback from the heated-windshield 
area is of little consequence, no effort is made to evaporate all 
the water intercepted by the windshield; such a requirement would 
impose an exorbitant heating load in some cases. 

Using the eqviations of reference 5 and the symbols of this 
report the rate of heat loss due to evaporation is 


^2 


0.622 h Lc 


°pP'o 


(8) 


Heat Loss Due to Warming the Impinging Water 


Eeference 5 gives the basic eqtiation for the dissipation of 
heat to the water that is intercepted by the windshield, but a 
kinetic heating term is added in this report since the loss of 
energy by the water droplets asstanes an appreciable value at 
airspeeds of. 300 miles per hoiur or greater. 


where 



and M is the average weight rate of water impingement per unit 
area of windshield and is given by 

M = 0.225 -3- Uom — (10) 

100 °^ A 


Since the Streamlines and drop trajectories for various types 
of windshields are not now available to calculate the efficiency of 
water impingement t),. data for spheres or flat ribbons can be 
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utilized to give approximate impingement values. Curves presenting 
efficiencies of water impingement for cylinders, spheres, and 
rihhons are given in reference 4 together with a discussion of 
their use. 


Heat Loss Due To Radiation 


The heat loss due to radiation from the windshield svirface to 
the surrounding atmosphere is usually of small magnitude. It can 
he calculated from the Stefan— Boltzman. equation: 



Summation of Heat Losses 

5 . . . ’ 

Combining equations (1) to (11), the complete equation for the 
dissipation of heat from a. windshield surface in conditions of 
icing may he written . 


= h t 3 -to ^.832 r 


2-1 


100) J 


+ 0.622 hL, 


(¥K) 


+ Me 


Pw 




+ 0.173 € 


I- V 4 . .4 

( ./ZL.) 1 

L\ 100/ Vioo/ J 


(12) 


For flight speeds under approximately 200 miles per hovir, 
equation (12) may he simplified hy neglecting the heat loss due to 
radiation and the kinetic heating of the impinging 'water. Regrouping 
the terms to segregate the primary heat losses, which are due to 
convective heat transfer and evaporation, from the secondary effects 
of heating the impinging water and kinetic heating of the air, equa- 
tion ( 12 ) becomes 


q - hX(tg-to) + M(ts-to) - 0.832 hr 



( 13 ) 
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vhere 

T - 1 + 0.622Lg f 63—00 ^ 

CpPo Ug-to/ 

A similar grouping was used in references 

DESCPIPTION OF APPAEATUB 

Hie research reported herein was conducted in flight using a 
C— 46 cargo airplane (reference 3) modified to incorporate an MCA 
thermal ice— prevention system which permitted continuous operation 
in natural icing conditions, The airplane is shown in figure 1 as 
equipped for the winter of 1945r-46, and in figure 2 as operated in 
the winter of 1946-47. 

Meteorological equipment was installed to measure the free— 
stream air temperature and the water drop size, drop-size 
distribution, and liquid water content in icing clouds . A shielded 
thermocouple connected to a millivoltmeter was used to measure the 
air temperature. Eotating cylinders of different diameters, 
mounted coaxially on a single shaft, were used to determine the 
liquid water content, drop size, and drop-rsize distributions. The 
maximum drop size was determined from the area of water impingement 
on a single, nonrotating cylinder. A detailed description of the 
meteorological instruments, their use, and typical meteorological 
data obtained are presented in reference 9* 


( 1 ^) 


5 and 8. 


Internally Heated Windshields 

Three general types of aircraft windshields were tested to 
provide fundamental heat— transfer data: flat— plate, flush, and 
V— type windshields. The first two types were tested during the 
1945^6 operations and the third during the 1946-47 operations. 

Each windshield— panel test section was heated by electrical power, 
and the installations are shown in figures 3> and 5. 

Flat— plate windshield .— This installation was intended to supply 
heat— transfer data for two-dimensional flow over a flat plate inclined 
at different angles to the air stream. It was considered that the 
data would be applicable to the design of flat— pane windshields such 
as those common in fighter and transport aircraft. The test panel 
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was located on the right side of the fuselage of the airplane as shown 
in figure 3* panel was mounted in such a manner that the angle 

■between its surface and the relative wind direction could be varied. 
Provision was made to fix the panel in one of three positions: 30 °, 

45°, and 6o°, measiired from the tangent to the fuselage at the hinge 
line, 

Ihe over-all surface dimensions of the flat plate were 13 by 
l8 inches, and of this area 13 by l6, 5 inches, or 1,49 sq^uare feet, 
were uniformly heated. Ihe heating elements were thin metal ribbons, 
and were capable of dissipating 2500 Btu per hour per square foot of 
panel sxirface when connected in a 24r-volt. circuit. For the base of 
the panel, a 1/2— inch sheet. Of fabric— laminated plastic was used. 

The heating ribbons were^CeIaented : to '' this base and then a 1/64— inch- 
thick sheet of laminated plastic was. cemented on top to present a 
smooth, electrically insulated ■ surf ace tjO the air stream. Centrally 
located in the base material was a -.4,5— inch-square heat meter 
(calibrated thermopile, reference io), l/64— inch— thick, to measure 
heat flow away from the test surface.;- Electrical power was supplied 
from an alternator .through ■a-120/24— volt transformer, and the 
voltage impressed on the test section was controlled by means of a 
carbon-pile rheostat. 

Temperatures on the suld'ace of the flat plate were measured at 
10 points by iron-constantan thermocouples which were 0.002. inch 
thick. The thermocouples were recessed in grooves 0,003 inch deep 
in the l/64-incb-plastic outer sheet, and their Junctions located as 
shown by figure 6. , The outer surface was thinly painted. to provide 
a smooth surface and to protect the thermocouples. Additional 
thermocouples were located on the back surface of the plate and at 
the heat meter. 

All temperatures were recorded by a self-balancing automatic- 
recording potentiometer with an estimated over— all accuracy of ± 3° F, 
This degree of accvtracy was made possible by the provision of several 
calibration thermocouples, with temperatvires knoim to ±1° F, as a 
check on the recording potentiometer.. The potential of the heat 
meter was determined manually with an additional potentiometer. By 
measuring the tempera tvire at the heat , meter, an accuracy of 
measiurement of heat losses to ±10 percent was possible. 

For the determination of the boundary— layer— velocity profile 
in clear air, conventional pressure rakes were utilized and the 
pressures were recorded by photographing a manometer board. 
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Ice formations on the panel were photographed with a 
16-mlllimeter cinematic camera which was remotely controlled. Ihe 
camera was- located within the fviselage about 3 feet forward of the 
flat— plate panel, and was fitted with a 90° angle lens which 
protruded into the blister shown in figure 3» 

Flush windshield .— This windshield was installed flush in the 
copilot's panel which, at that time (1945—46), was flush with the 
fuselage contovirs. (See figs. 3 and 4.) The test surface was 
almost sqmre, though slightly curved, and covered an area of 
1.16 square feet. Construction of the panel consisted of three 
sheets of glass bonded together with sheets- of plastic. Heating 
elements of very fine resistance wires were imbedded between the 
outer glass layer and the plastic inner layer. Ihese wires were 
0.0012 inch in diameter and spaced 0.6 inch apart. Bie panel was a 
ccamnercial product and was designed for a maximum heat dissipation 
of 1000 Btu per hour per square foot. Power was supplied by direct- 
current generators and the impressed voltage was varied with a 
rheostat. Thermocouples 0.002— inch— thick were cemented to the inner 
and outer surfaces of the test area and covered with clear spar 
varnish. The locations of the thermocouple junctions are shown in 
figure 6. 

Ice formations were photographed with a box camera having flash 
attachments and located in the pilot's compartment. Temperatures 
and boundary— layer-rake pressures were recorded by the same instru- 
ments used for the - flat— plate windshield. . 

• V— type- windshield .— For the 1946—47 operations the flush wind- 
shields were replaced by the V— type windshield shown in figure 5» 

The windshield center post was sloped at an angle of 57° (measured 
from the horizontal during cruising flight) and the included angle 
of the V was 86° measured in the horizontal plane. The fabrication 
of the pilot's and copilot's panels was identical. and consisted of 
three sheets of glass bonded together with interposed sheets of 
plastic (fig. 7). Heating of the panels was attained through a 
transparent, electrically conductive surface coated on the outer 
sheet of glass on the surface adjacent to the plastic inner layer. 

The panels were a commercial product and were capable of dissipating 
2000 Btu per hour per square, foot. Althoiogh the panels were 
trapezoidal in shape, the heated area was rectangular with the bus 
bars applied. the full length of the rectangle along the top and 
bottom edges to insure uniformity of heat distribution (fig. 7)« T^ie 
heated area of each panel was 14.5 by 25 inches or 2.52 square feet^ 
Power was sppplied by an alternator and the voltage was controlled 
with a variftbjLe transformer. Thermocouples 0.002— inch— thick were 
cemented to t^e inner and outer surfaces of both windshields and 
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covered with clear spar vanish. Locations of the thermocouple 
junctions are shown in figure 6. 

Photographs of ice accretions on the windshields were taken 
with a camera having a synchronized flash attachment. Temperatures 
were recorded in the same manner as descri'bed for the flat— plate 
windshield. 


External Air-Heated Windshields 

Flush windshield ,^ The original windshield heating system 
installed in the C— 46 airplane included hoth an internal douhle— 
panel heating system and an external— discharge heated— air heating 
system as discussed in ref erences 3 Q-nd 11. The original system 
was modified as sho^m in figure 8 for the 1945-46 operations. The 
inner heating system was removed, the heated— air discharge slot for 
the outer system was decreased to three— sixteenths inch, and the 
extent of the discharge' slot was reduced to. a width of 1 foot at 
the hottom of the pilot *s windshield, starting at the center post. 

The heated air was supplied hy exhaust-gas heat exchangers, and the 
air flow rate was measured with a 4— inch venturi meter. A butterfly 
valve was installed to control the quantity of heated air flowing 
from the discharge slot. 

Windshield surface temperatures were measiared with thin iron— 
constantan thermocouples cemented to the glass and covered with 
several coats of spar varnish. Nine thermocouples were located on 
the outside surface and three were located on the inside surface, 
as shown in figure 6, . Thermocouples were located at the venturi 
meter and just inside the discharge slot to measure the heated— air 
temperatures. Two temperature rakes were installed on the windshield 
outer surface for the determination of the temperature profile in . 
the heated-air jet in a few clear-air flights. These rakes consisted 
of nine thermocouples at different heights above the windshield 
surface up to 1-1/2 inches. All thermocouples were connected to a 
recording potentiometer,* and the estimated over-all accuracy was 
about ±3^ F. 

A pressure rake was installed for the determination of the 
heated-air Jet— velocity profiles in clear-air flights. The rake 
consisted of l8 pressure probes mounted in the form of a triangular 
prism with its base on the windshield surface. Three of the 
pressure probes were static tubes and the remainder total— pressure 
•tubes. These pressures, along with the venturi meter pressures, 
were recorded by photographing an alcohol manometer board. 
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IVo methods vere employed to obtain an indication of the 
qtaantity of heat flow from the external jet through the windshield 
into the airplane cabin. The first method consisted of the calcula- 
tion of the heat flow based on the measured temperature difference 
between the windshield surfaces and the known thermal conductivity 
of the windshield materials. In- addition to this calculation, a 
k. 5-inch-square heat meter was located on the inner surface as shown 
by figure 8. The potential across the heat meter was read by an 
indicating potentiometer, Ihe estimated accuracy of the. heat flow 
as determined with the heat-meter installation is approximately- 
±15 percent. 

Y— type windshield .— For the 1946-47 operations with the V— type 
windshield, one 3/18— inch external discharge slot was located along 
the bottom edge of the pilot's windshield and another along the 
center post between the windshields with the heated air directed 
over the pilot's windshield. Each slot was 1 foot in length. The 
4-inch venturi meter was replaced by two 3^inch ventiuri meters 
equipped with thermocouples to measure the flow rate independently 
to the two separate discharge slots over the pilot’s windshield. 

The air— flow rate to each slot could be adjusted independently 
although a certain amotint of interdependence existed,- since both 
slots received air from the same source of supply. A booster blower 
was installed in the main supply duct to increase the total qviantity 
of heated air available to the slots . 

The fabrication and location of thermocouples installed on the 
pilot's panel were described in the preceding description of the 
V—type windshields. In addition to the thermocouples on the wind- 
shield surface, four thermocouples were installed just inside the 
center— post discharge slot, and six thermocouples were located just 
inside the slot at the bottom of the pilot's windshield. These 
thermocouples were located so as to measure the heated— air 
temperatures just before the air was discharged from the slots. A 
self-balancing potentiometer was used to record the tempera tiores . 

The pressure differentials at the venturi meters were indicated 
by standard differential pressure gages and recorded by a standard 
MCA 60-cell recording manometer. 


TEST PROCEDURE . 

The clear— air flights were conducted at the Ames Aeronautical 
Laboratory, Moffett Field, Calif. In the case of the icing flights, 
the Laboratory -served as a main base and operations in search of 
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lolng conditions vere conducted over a considerable area of the 
United States, The 1945-46 operations, consisting of over 40 
flights, were conducted mainly in the Pacific Northwest and North 
Central States, The 1946-47 operations included this same area 
and, in addition, some icing conditions were encountered in the 
Central Southern states. 

The flights in clear air were made to check or modify the 
theoretical equations for convective heat-transfer coefficients 
and to measwe the degree of kinetic heating experienced by the 
windshields. 


Icing :Flights * 

The test proced\ires employed were substantially the same for 
each test windshield. Therefore, the following discussion applies 
to the operations of both winters. When icing conditions were 
anticipated, the windshield power was txumed on or, in the case 
of the heated-air jet, the valves were opened. Before entering 
the icing condition, the heat to the windshields was adjusted to a 
setting greater. than the expected requirements. After a stabiliza- 
tion period of several minutes the recording instruments were started 
and the cloud was entered. After the completion of one run, if no 
ice formed, the heat supply would be decreased and the equipment 
allowed to stabilize for 10 minutes before the next set of data was 
taken. The heat input to each windshield was reduced in this manner 
to approximately the minimum amount required to keep the heated 
surface free of ice accretions. For the pilot’s and copilot’s wind- 
shields, the value of minimum heat required was established by 
visml observation of the start of ice accretions on the windshields. 
Since the swface of the flat-plate panel could not be observed by 
the windshield engineer, the outer— surface temperatures were 
regTilated by changing the heat input so that the lowest temperature 
wovild be just above freezing. 

The pilots always made an attempt to hold the altitude and 
airspeed constant dtiring a specific test rvin, but severe turbulence 
sometimes made this impossible. The dwation of a test run varied 
from a few minutes in the top bf a cumulus cloud to several hours 
in an e^rtensive stratus layer. For a few runs in the smaller cloud 
formations, the airplane was flown in a circle in order to remain in 
the cloud long enough to obtain a complete set of meteorological and 
heat-transfer data. None of this latter type of data has been 
presented \inless the conditions were fairly stable for at least 
5 minutes after an equipment-heating stabilization period of 10 
minutes. 
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In addition to the fli^t and meteorological conditions, a 
continuous record vas kept of the following windshield data: heat 

supplied, inner- and outer— surface temperatures, photographs of ice 
formations, and the heat losses inward as indicated by heat meters. 


Clear— Air Tests 

Clear-air flights were conducted in level flight at various 
altitudes and airspeeds to cover the conditions encountered during 
the icing operations, Ihe test flights were made in a direction 
away from the sun to minimize the effects of solar radiation. 

First, flight data were taken without any power supplied to 
the windshields. By measuring the surface— temperature rise above 
free— stream temperature a basis was established for comparison 
between the indicated recovery factor of equation (4) , and the 
theoretical value Pr for turbulent flow. 

A check for possible significant edge losses was made by varying 
the heat supplied to the windshield panels while flying at a constant 
airspeed and. altitude. , Any variance in the calculated external heat- 
transfer coefficient after correction for inwajrd heat losses was an 
indication of the edge losses, since the effect of changes of 
tanperature in the boimdary layer on the heat— transfer coefficient 
can be considered negligible for the temperatiire ranges concerned. 

For several airspeeds and altitudes, heat was- supplied to the 
windshield panels while the temperatures, heat losses, and velocities 
were recorded. These data were necessary for the experimental 
determination of the convective heat-transfer coefficient h, and 
the establishment of effective values of s or I (eqxiations (6) 
and (7)). 

During the clear-air flights after the 1945-46 operations, a 
survey of the boTindary layer at the center of each test windshield 
was made with a pressure rake. The data were used to determine the 
boundary— layer thickness and the variance of the local velocity 
outside the boundary layer with changes in airplane velocity. Both 
preqqure ajicL temperature boundary— layer surveys were made in the 

heated^ir jet over th© flush windshield installation. 

% - / 
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■RESULTS OF TESTS 
Internally Heated Windshields 

Local air yelocities and convective heat^transfer coefficients ^ 

The values of local* air velocity just outside the boundary layer at 
the center of the flat— plate panel and the flush windshield are 
presented in figure 9 for various airplane speeds and altitudes. In 
figures 10, 11, and 12 the convective heat-transfer coefficients, 
based on the measured surface temperature at the center of the flatp- 
plate and the known heat flow, are presented* The coefficients for the 
flush panel are not presented because of difficulties experienced 
with the measurement of the surface temperature . In the case of the 
V— type windshield only . a ‘few values of the heat— transfer coefficient 
obtained in clear air are available. These data were secured 
between the icing encoixnters and, since the data were recorded at 
different flight altitudes, they do noi form a satisfactory basis 
for cui*ve plotting in the same mariner as figures 10, 11, and 12. 

This, information is presented in the following table: 


Ifest 

no. 

Pressure. 

altitude, 

(ft) 

True 

airspeed 

(mph) 

1 

Free-air 

temperature 

(°F) 

Windshield 

aur£cco 

tfemperature 

(Op) 

h 

Btu/hr, 
sq. ft, 

(%) 

1 

■ 12,000 

178 

20 

73 

23. U 

2 

10,800 

200 

18 

65 

28.5 

3 

9,800 

17h 

39 

73 

23.0 

k 

. 7,800 

175 

-30 

66 

26.5 

. 5 

6,000 

170 

35 

70 

25.2 


Meteorological and heat-transfer data during icing conditions .- 

Thbles I and II present .the meteorological data for the two winter *s 
operations which have been selected for discussion in this report. 
Thesri data represent only a small portion of the meteorological 
information recorded, but have been selected as the tests which 
supply the most satisfactory combination of simultaneous meteoro- 
logical and windshield heat-txansfer data. A complete tabulation 
of all the meteorological data obtained in the 1945—46 operations 
is presented in reference 9 , which also includes a detailed 
discussion of the procedures employed to obtain the test restilts.. 
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!I3ie rotating-cylinder method utilized to measure the liquid water 
content and meanr-effective drop-size values, presented in tables I 
and II, provided a 1-minute average of the icing conditions, since 
this was the time of exposure of the xjylinders.^ The windshield 
heat- transfer data required about ' 2 minutes to record and the 
rotating cylinders were exposed at some time during that interval, 
Ihe meteorological data of icing conditions k, 5, and 6 of table I 
were obtained prior to the Installation of the rotating cylinders, 
hence the maxi mum drop-size range encoimtered during the heat- 
transfer test interval is presented. 

Ihe heat— transfer data for all three of the internally heated 
windshields in conditions of icing are presented in table IH, All 
the heat flow. values presented in the table do not represent the. 
minimum requirement for the particular icing condition, since cases 
for which satisfactory meteorological and heat— transfer data were 
obtained together did not always occur at the minimimi value of heat 
flow. The surface temperatures presented in table HI are the 
arithmetical average of the recorded values for the several thermo- 
couples on the surfaces . Photographs and sketches showing the ice 
accretions on the heated windshields and surro^lnding surfaces for 
several of the conditions encountered are presented in figure I3 
to 23. 


External Discharge Windshields 

Temperature and velocity gradients .— The temperatures of the 
windshield surface at various distances from the heated-air discharge 
slot for the flush windshield installation are shown in figure 2 k, 

The temperature and velocity gradients in the Jet for the same 
windshield installation are presented in figures 25> and 26, 

Icing tests .— Table IV presents the windshield (flush) surface 
temperatiires measured during operation of .the external discharge 
jet in icing conditions. An indication of the ice— prevention action 
of the external discharge jet for the flush windshield installation 
can be seen in figure I5. The effect of the V— type windshield 
external jet (pilot's side df windshield) is shown in figure 20. The 
incomplete removal of ice accretions shown in this figure was typical 


^Mean^ffective drop size for a cloud sample is defined as the 
diameter of a water drop for which the amount of water existing iri 
water drops larger than that drop equals the amount of water in drops 
smaller than the drop. 
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of the operation of the external discharge system on the V— type 
windshield, ISie temperatxire limitations of the windshield panel 
precluded the use of additional heated air in an attempt to remove 
the residual accretions., 

DISCUSSION 

Prediction of Heating Requirement for 
Internally Heated Windshields 

Since the main ohjectlveof this report is to establish 
fundamental design equations for the calculation of windshield heat 
requirements in specified icing conditions,, equation (I3) will he 
discussed in the light of; the test data’ to establish its usefulness 
for this purpose. Ihe ap^plicability of this equation for the 
determination of heating requlranents can best be evaluated by a 
comparison of calculated and measuired values of heat flow from the 
windshield surface, for the same surface temperature, and for 
specified flight and icing conditions, Before such calculations 
can be made, however, means -must be established for the determina- 
tion of two components of eqmtion (I3) ; namely, the convective 
heat— transf er coefficient h • and the weight— rate of water impingement 
M. Consequently the following discussion considers first the 
utilization of the flight— test data to establish means for the 
evaluation of h and M with a degree of accuracy suitable for 
design purposes. Althotigh the methods discussed are directly 
applicable to the test windshields of this report, an effort was 
made to generalize the technique in order to provide a design 
procedure applicable to other configurations and icing flight 
conditions. 


Evaluation.,of ■ the, convective heat-transfer coefficient .- In the 

determination of the value of the convective heat— transfer Coeffi- 
cient- h, the pertinent question' is the suitability of equation (6) 
for the calcvilation of the coefficient. A comparison of measured 
and computed values of the convective heat— transfer coefficient h 
for the outer surface of the three internally heated test wind- 
shields therefore will be made. 

Considering first the evaluation of h for flat— plate panels, 
figures 27, 28, and 29 present a comparison of the measured values 
of h at each panel angle and at various altitudes and airspeeds, 
with values calculated by equation (6) . The curves representing the 
measured values of h have been teken from figures 10, 11, and 12. 
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Equation (6) is used for comparison rather than equation (7), since 
it has heen derived for the purpose of calculating the coefficient 
at a given point on the surface and the test data apply to the center 
of each panel. The values of Ui used in equation (6) for the 
calculation of the coefficients presented in figures 27, 28, and 29 
were the local velocities at the center of the panel as given in 
figure 9. TWo values of s were used in the calculations Si 
representing the distance from the panel hinge to the center of 
the panel, and sg the distance from the center of the panel to 
the stagnation point at the nose of the airplane fuselage. 

In the case of the V— type windshield, the .calculated values 
of h from equation (6) ard compared with the test values (see 
table presented in Results) in the following table: 


Item 

no. 

h 

measiired 

h calculated 

(equation (6)) 

Bi = 0.67 ft ! 

I 82 = 5.0. ft 

1 

21.3 

29.1 

19.5- 

2 

25.9 

33^2 

22.1 

3 

20.8 

29.6 

19.8 

k 

24.0 

32.4 

21.6 

5 

22.7 

33.1 

22.1 


In the computation- of h by equation (6) the value of the local 
velocity over the windshield was not known and hence the airplane 
velocity was used. 

■ An examination of figures 27, 28, and 29 and the table Just 
presented indicates that the use of equation (6) provides reasonably 
accTirate values of the convective heat— transfer coefficient at the 
center of the test windshields. In terming the accuracy of eq\ia— 
tipn (6) as "reasonable," consideration has been given to the fact 
that no data other than the length a.nd location of the panel and the 
local velocity were employed and the agreement is considered reason- 
able for such an approximation, A more exact determination of the 
value, of h would require the application of somewhat lengthy computa^ 
tion methods to measured values of the boundary— layer profile, or 
the installation of a heating plate in the windshield surface and the 
procurement of actual test data. The disadvantages of these means 
must be weighed against the inaccuracy of the general application of 
eqmtion (6) in determining the method that will be employed to 
evaluate h in an.y futtire design computations. The agreement 
shown in figures 27,- 28, and 29 between the measured values of h 
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and the veilues con 5 )uted hy equation (6) is about the same, 
irrespective of the value of s employed. In most instances 
the two ccanputed curves bracket the experimental curve, with the 
lower value of s providing the larger values of h. For design 
purposes the use of the lower value of s should provide heat- 
transfer coefficients 'which wotild probably be somewhat larger 
than the true case, but would be conservative. 

Ihe spacing of the heating wires in the flush windshield 
(0.6 in. between wires) resulted in considerable nonuniformity of 
the surface "tenperature as can be observed in the ice— removal photo- 
graphs (fig. 15 ) • As a result .the atteupts to, calculate the 
convective heat— transfer coefficient for the surface based on the 
measured heat flow and' surface- temperature, rise produced question- 
able £ind erratic values. Consequently a comparison of measured and 
cQEputed values of h for the flush windshield has not been presented. 

In an atteupt to establish the accuracy of equation (I 3 ) for 
the prediction of the heat requirement, it is desirable to utilize 
the most accurate available values of the various components of the 
equation. Thus although the foregoing discussion indicates that 
equation (6) provides reasonably accurate values of h (which can 
be used when more precise methods are not available or are 
considered too conplicated for the degree of accuracy desired), a 
more acctirate determiination of h would be desirable for use while 
checking the general equation. One possible solution would be to 
utilize the teat values obtained in clear air, but this system has 
the disadvantage of inflexibility since corresponding clear— air 
data were not available for all the icing tests. The alternative 
euployed was to determine the value of s in equation (6) which 
would provide the best average agreement between calculated and 
measured values of h for the test range of airspeeds and altitudes. 

13ie determination of an average s for the case of the flat- 
plate panel at an angle of 45*^ will be used as an example of the 
procedure followed. In figure 11, 12 test points are plotted as 
measured values of h. The data corresponding to each test point 
(h, T, U, and 7 ) were inserted in equation ( 6 ) and the value of s 
computed for each point. The arithnetical average of the 12 values 
of s was calculated to be 1.68 feet. For the determination of the 
value of h at the center of the fla 1 >-plate panel when set at 45° 
for different velocities and altitudes, equation ( 6 ) becomes: 

h = 0.51 Tav°'® ( 1 ^ 68 )°’° = 0.46 Tav°-® (Ui7)°'® 


(15) 
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This same system vas eu^loyed for the other angles of the flat- 
plate panel, and for the V— type windshield, with the following 
results: 

Panel angle = 30°, s = O.37 foot, and h = 0.62 Tav°’® (16) 

Panel angle = 60°, a = I.93 feet, and h = 0.45 (Ui7)°'® (I7) 

V— type windshield: 

s =1.26 feet, and h = 0.4l Tav°*® (Uo7)°*® (I8) 


Evaluation of the rate of water impingement .— Assignment of a ; 
value to the factor • M of eqviation (ij.3) for specific design purposes 
is a problem reading which very little information is available. 

The most recent and extensive treatment of water-drop trajectories 
around several generalized objects is presented in reference 4. This 
investigation was undertaken primarily to provide a method for 
estimating the rate of water impingement on an airfoil and, there- 
fore, most of the report concerns the interception of water drops 
by. cylinders which are assumed to represent the forward section of 
the various airfoils. Some trajectory calculations, however, are 
presented for the cases of a ribbon (flat plate normal to the 
direction of air flpw) and a sphere. These calculations are used 
in this report to predict the rate of water iag)ingement on the flat- 
plate panel (ribbon) and the flush windshield (sphere). 

Considering first the case of the flat— plate panel, the 
assumption was made that the rate of water impingement on the panel 
woxild be equal to the rate of impingement on the projected area of 
the panel considered as one— half of a ribbon. (See fig. 30. ) The 
efficiency of water impingement q and the weight rate of water 
iDq)ingeinent w for various drop sizes have been computed for the 
three panel angles at one flight condition and are presented in 
table V. It is of interest to note that for drop diameters greater 
than 30 microns, as the panel angle is reduced (panel becoming more 
flush with the fuselage) the impingement efficiency is increased. but 
the wei^t rate of water interception is decreased. Since this 
latter quantity is the item of greater significance from a heat- 
requirement standpoint, the desirability of low panel angles is 
evident. Optimum benefit of this advantageous feature would be 
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expected to occur with flush windshield installations. The flush 
windshield in the C— 46 airplane dxiring the 1945 t- 46 winter provided 
confirmation in this respect, since windshield icing was encountered 
in only 8 out of 30 actual flights in icing clouds. Figure 16 shows 
the lack of ice accretions -on the flush windshield after a flight in 
a cloud which produced the accretions shown on 3/8-inch rods located 
Just helow the pilot's. and copilot’s windshields. 

A second item of interest which can he noted in table V is the 
fact that the maximum calculated efficiency of water impingement is 
80 percent,' and this value corresponds to exceedingly large drops 
for an icing condition. This restilt is somewhat surprising, and 
possibly in error, since it wovild appear that the presence of the 
fuselage forward of the windshield would cause a concentration of 
drops near the surface and produce impingement efficiency values of 
large magnitude, possibly greater than 100 percent. Some attempts 
were made to obtain an indication of the, rate of water impingement 
on 'the C— 46 windshields durihg the opera -bions, but nothing satis- 
factory resulted, 

- To obtain an estimation of M the rate of water impingement 
on the flush windshield for specific design conditions, the 
assumption was made that the impingement would be the same as that 
on a portion of a sphere having a diameter equal to the fuselage 
maximum diameter. (See fig. 31.) The manner in which' the water- 
drop trajectory calculation method for a sphere presented in 
reference 4 was applied to the configuration shown in fig;u?e 3I is 
discussed in- detail in Appendix A of this report. The calculated 
rate of water impingement on the flush test panel, based on the 
procedure discussed in Appendix A, is presented in table VI for 
various drop-size diameters.. The same flight conditions employed 
for the flat— plate water interception calculations (table V) were 
again assumed to apply. 

An interesting conclusion resulting from an inspection of 
•table VT is that, for the conditions presented, drop diameters of 
50 microns are required in order to achieve a value of % of 
sufficient magnitude to reach the test panel. This conclusion is 
in agreement with the observations of water impingement on the 
C— 46 flush windshield in icing conditions. The drops were observed 
to strike the windshield only d\u*ing flights -for which the 
meteorological data indicated .the presence of mean effective drop 
size of 30 microns in diameter, or larger. The fact that the rate 
of water impingement on the flush panel as presented in 'table 71 
decreases with increased drop size is the result of the assumption 
(Appendix A) that •the water is uniformly distributed over the area 
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of iiqpingeinont ^ and means that the -eirea of impingement increases at 
a more rapid rate than the weight rate of water iH 5 )ingement.. This 
possible error in the value of the wei^t rate of water in^jingement 
is considered to be of negligible magnitude in view of the amftn 
amount of water involved for the drop— size range assumed. 

Establishment of applicability of eq,mtion (I 3 ) for the 
prediction of windshield heating requirements .— The applicability 
of equation (I 3 ) will now be demonstrated by a coD5)arison of calculated 
and measrired values of heat flow in conditions of icing for the flat- 
plate, flush, and V-type windshields. The assuB 5 )tions and method of 
calculation followed are presented in Appendix B. To illustrate the 
procedure employed, sai^le calculations are presented in detail for 
the flat plate and flush windshield for icing condition 6 , tables I 
and III, An illustrative calculation for the V— type windshield is 
not included, since the procedure followed was identical to that 
presented for the flat-plate panel. 

Ihe calctilated veilues of heat flow are compared with the measured 
values for the flat— plate _and V— typo installations in tables VII and 
VIII, respectively. The fact that the calculated values of q are 
lower than the measured values in table VII and higher in table VIII 
may be attributed, in part, to the different velocity basis used for 
the determination of h in each case . The calculated values in 
table VII should be more nearly correct since the velocity enq>loyed 
was the local velocity over the panel, while in the case of the V— type 
windshield (table VIII) the airpl^e speed was used since the local 
velocity over the windshield was not known. 

The few instances in which icing of the flush test panel 
occurred provided very limited data for a conq)arison of the calculated 
and measured values of the heat flow. The one case presented as an 
illustration in Appendix B represents the only reliable data, and 
the calculated value was percent greater than the experimental 
veilue. 


One possible reason for the large disagreement, in some cases, 
between the measured and calculated values of heat flow may be 
traceable to the assunption, used for all of the calculations, that 
all the drops were of the same size. The calculation of the heat 
flow from the flat-plate panel for icing condition 12, table VII, 
will be used to illustrate the . inf luence of drop-size distribution on 
the heat flow. The calculated veilue of 89 O Btu per hour, square 
foot presented in table VII, icing condition 12, was based on a 
\miform drop size of 19 microns. The rotating-cylinder data, however, 
indicated a drop-size distribution defined as type E in reference 4. 
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Tfe.'ble IX presents the resxilts of calculating the weight-rate of water 
impingement on the flat-rplate panel, assuming the water in the free 
stream exists in drop sizes corresponding to distribution E, On this 
basis the original weight rate of water impingement (1.91 lb per hr) 
is increased to 8.l4 pounds per hour. . Although this is a sizeable 
increase in the amount of water intercepted, only q3 (equation (9) ) 
of equation (13) is affected by the increase, and the revised value 
of q becomes 1035 Btu per hour, square foot, which, in this case, 
provides better agreement with the measured value. 

•- The -assumption of miform drop-size distribution in the calcula- 
tion of q ^y also be the cause of the large disagreement between 
the measured and calculated values of q in Icing conditions 1 and 
7, table VII. Based on the assumption of no drops present larger 
than the mean effective diameter, the weight— rate of water impinge- 
ment in these two cases is foimd to.be zero.' The value of q in 
table VII, therefore, has been derived for these two conditions on 
the assumption of a clear panel with no evaporation occurring on 
the surface'. If drops of a diameter large enough to strike the panel 
had actually been present in sufficient quantity to wet the panel, 
even partially, the calculated value of q would be considerably 
increased and better agreement with the measured value would result. 
Although these discussions are by no means conclusive, they do 
indicate the importance of drop— size distribution and maximum drop, 
size. This factor shovild be given consideration in future design 
calculations. 

Based upon an examination of the data in tables VII and Vlli, 
the . calculation for the flush windshield heat flow in Appendix B, 
and a consideration of the various influencing factors which cause 
disagreement between the calculated and measured values of heat flow, 
it Is concluded that equation (13) will provide predicted values of 
heat flow with an average accuracy of, 15 percent if the entire surface 
is assumed or known to be completely vetted and the local velocity 
over the windshield panel is employed in the calculations. The 
accTAracy of the eqmtioA is greatly dependent upon the accxxracy of 
the values of convective heat— transfer coefficient employed. Tho 
determination of experimental values of this quantity for the 
particular windshield configuration proposed would be of considerable 
aid in the accurate prediction of the required heat flow, although 
in the absence of such data equations (6) and (T) will supply a 
reasonable approximation. The effect on the accuracy of eqxiation (13) 
produced by an error in the estimation or calculation of the weight 
rate of water impingement is not large, provided sufficient water is 
intercepted by the windshield to completely wet the surface. Bates 
of water impingement less than the value required to wet the surface 
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affect the value of qa (equation ( 8 )), since evaporation over the 
entire axirface is not realized. For rates of impingement above the 
quantity required for completely wetting the svirface, the value of 
qa remains constant and the only term of equation (I 3 ) affected is 
the value of q 3 (equation ( 9 )) which is an item of secondary 
magnitude. 

General considerations of windshield heating requirements .— 
Having shown that equation (I 3 ) is, applicable for the calcvdation of 
windshield heating requirements, provided the various components of 
the equation can be evaluated, it is of interest to utilize the 
equation to investigate the heat requirements for- various windshield 
configurations in icing conditions which might be selected as design 
requirements. An indication of the meteorological conditions , 
corresponding to typical, or average severity, icing and to the most 
probable maximum icing to be expected in all weather transport , 
operations was obtained from a -review of reference 9» For each of 
the two general cloud types (cumulus and stratus) a liquid water 
content was selected which corresponded to a moderate icing condi- 
tion. In the case of the stratus cloud, the value assumed was 
,0.5 gram per cubic meter, while the corresponding value for the 
cumulus cloud was 1,0 gram per cubic, meter. Based on the da'ta in 
figures 6 and 9 of reference 9^ values of 15® F, free— air temperature 
and 15 microns, drop diameter, were selected for the stratus cloud. 
The corresponding values for the cumulus cloud were 0° F and 20 
microns. , ' . 

The heating requirement for the flat— plate, flush, and Y— type 
windshields in the above icing conditions was calculated by the 
method presented in Appendix B for airplane speeds of 150 and 
300 miles per hovir, and the results are presented. in table X> From 
this table the previous empirical heating requirement of 1000 Btu 
per hour per square foot of windshield surface (reference 1) is 
seen to provide adequate protection for V— type windshields in 
moderate stratus for flight speeds as great as 3 OO miles per hour. 

In the case, of the moderate cumulus cloud, however, this quantity 
of heat flow would not be' adequate even- at a speed of I 50 miles an 
hour. In transport aircraft operations, the most critical need for 
windshield ice protection occurs during flight in the congested area 
surro'unding the air terminals. Any icing condition encountered would 
be of the stratus type which would tend to reduce the heating 
requirement, • For transport operation in mountainous areas, where 
cumulus and stratocumxilus; predominate and high cruising altitudes . 

. (hence low free-air temperatures) are the rule, a heat require- 
ment of from 2000 to 2500 Btu per hour per square foot, for speeds 
up to 300 miles per hour, is Indicated if complete and continuous 
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protection is desired. Certainly complete visibility at all times 
is desira'ble for military aircraft. It is of interest to note that 
a heat supply of iCi^^Btu .per hour per square foot appears adequate 
for a flU'Sh windshield in all of the icing and flight conditions 
assumed in table X, 

In considering the individual contribution of each of the 
various meteorological factors ■v^ich constitute icing conditions 
to the total heat requirement , the powerful effect of changes in 
the value of free-air temperature 1s noted. This dominating effect 
is traceable to the fact that the surface evaporation is greatly 
influenced by temperature 'changes, and this te:m qs of the general 
equation is one of the major components. As an illustration of this 
effect, consider the increase in' heating requirement (table X) for 
the flat-opiate windshield from 5^0 to 1040 Btu per hotir> square 
foot as a result of (l) increasing the drop diameter by 5 microns^ 

( 2 ) doubling the liquid water content, and ( 3 ) decreasing the free- 
air temperature by 15° F. \Jhen equation (I 3 ) is utilized, it can 
be shown that the 500 Btu per hour per square foot increase is 
composed of a 13 — percent increase due to increased drop size and 
liquid water content, and 87 percent due to the change in the free- 
air temperature. This fact leads to the conclusion that once the 
■windshield heat requirement has been established for a specified 
icing condition, and when the condition is assumed to be' severe 
enough to completely wet the windshield surface, changes in drop 
size and liquid water content will not change the heating require- 
ment apnreciably, A change in free— air temperature, however, will 
have a very noticeable effect on the heating requirement. It should 
be noted that this conclusion is at variance with wing heating 
requirements because, in the case of wings, runback of the impinging 
water is not desirable and, therefore, sufficient heat is supplied 
to evaporate all of the water, wherever this is pmcticable. An 
indication of the amount of heat that would be required if this 
design procedure were applied to windshields can be obtained with 
equation (13) . Calculations with this.; equation, based on icing 
condition 6, table I-, -indicate that 'evaporation of all the water 
impinging on the flat— plate panel in this condition would require 
a heat supply of 19,400 Btu per hour, square foot. 


External Discharge Windshields 

Xr. the introduction of this report the statement irae made that 
measuT-ements Were taken in the jet of the heated— air external— 
dischar-.ge system in the hope that a basis for the establishment of 
dpplrical design, equations would result. Unfortunately, a review of 
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the velocity and temperatiire profiles measured in the Jet, of which 
figures 25 and 26 are typical examples, did not reveal any “basis for 
a rational analysis of the mixing of the Jet with the “boundary— layer 
air, and the prediction of the resultant surf ac ^temperature rise. 

A few items of interest, however, were noted during the investi- 
gation, One of these items was the large amount of heat supply 
required for ice protection by the external discharge method in 
comparison with internal heating of the surface,. For example, 
figure i5 illustrates a flight in which the windshield area under 
the discharge Jet was maintained clear with a heat supply of 20^000 
Btu per hour, which is an approximte imit flow of 10,000 Btu per 
ho\ir per square foot of cleared surface. In the same icing condi- 
tion, ice accretions were removed from the flush test panel with a 
heat flow of 5^5 Btu per hour per square; foot, (See fig, 15 and 
table in.) Thus the external discharge system required a heat 
supply approximately 20 times tMt required for the internally heated’ 
system, for the same degree of proteetioh in the same icing condition. 

The thermal inefficiency of the external-discharge system is 
apparently the result of rapid mixing of the discharge Jet with the - 
cold boundary-flayer air, with a resultant rapid decrease in' the Jet ' 
and windshield sixrface temperature as shown in figure 25. The 
surface tempera tizre data in this figure indicate a decrease from- 
l6o° F at the discharge slot to a value of 42° F at a distance of 
only 6 inches from the point of discharge. 

Satisfactory operation of the external— discharge system for 
the V— type windshield was not, obtained because of failure of the . Jet 
to flow across the entire surface. Figure 20 illustrates a typical 
ice— removal attempt, A small area along the bottom of the pilot’s 
windshield, and also at the center post, has been cleared, by' the Jet 
but the remainder of the panel is covered by an ice formation. 

Attempts to. raise the surface t^perature by increasing the flow 
rate and temperature of the heated air were limited by the temperature 
restrictions of the vinyl plastic at the .lower edge of the windshield 
(region of maximum temperature). The external— discharge system of 
windshield ice prevention appears to be a desirable installation only 
in those instances where (1) internal heating is not possible, 

(2) the discharged air will flow over the windshield without requiring 
the additional weight penalty of blowers, and (3) a large supply of 
heated air is available in the region of the windshield to be 
protected. 
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CONCLUSIONS 

Hie following conclusions are based on the analytical studies 
and test data of this report and should be applicable to windshield 
configurations and icing conditions similar to those investigated: 

1, Ihe coefficient of convective heat transfer for' the external 
surface of, a Y— type or flush windshield can be approximated with 
accuracy suitable for design purposes by the use of the established 
equations for turbulent flow on a flat plate. 

2, .The heat requirement for ice prevention on a flush or 

Y— type airplane windshield du:rtng flight in specified icing condi- 
tions can be calculated to ah>ccuracy 'of ; 20 percent. 

3, The complete and continuous prevention of ice accretions on 
the surface of a Y— type airplane windshield, for flight in moderate 
cumulous icing conditions at speeds up to 300 miles per hour, will 
require a heat flew from the s\irf ace of from 2000 to 2500 Btu per 
hour per square foot of surface. In the case of continuous flight 
for the same speed range in modeia.te stratus conditions, a. heat 
flow of 1000 Btu per hour per square foot should prove adequate, 

4, The complete and continuous prevention of ice accretions 
on the surface of a flush— type airplane windshield located well aft 
of the fuselage stagnation region, for. a speed range up to 300 miles 
an hour in stratus and moderate cumulous conditions, can be obtained 
with a heat flow of 1000 Btu per hour per square foot of siirface. 

5, The tendency of ice to. accrete on windshields which are 
installed flush with the fuselage contours is cpnsideiably less than 
that for Y— type windshields. 

6i The external discharge system of windshield thermal. ice 
prevention is thermally inefficient and requires a heat supply 
approximately 20 times that required for an internal system having 
the same performance. 

7. Windshield installations which conform to the fuselage 
contours are more adaptable to the use of the external discharge 
system than Y— t3^e installations because the heated jet will flew 
naturally over the windshield surface. 

Ames Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 

Moffett Field, Calif. 
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APPENDIX A 

Method Used for the Determination of the Rate 
of Water Impingement for the Plush Test Panel 

The cross sections of the C— 46 .airplane fuselage consisted of 
two circular sections separated by the floor line, the upper section 
being the greater. The entire flush windshield assembly was included 
between two stations which were located 23.5 and 50.5 Inches, respec- 
tively, from the nose of the fuselage and are designated as stations 
A and B, respectively, (fig. 3 I). 'Hie flush test panel was located 
in the copilot's windshield as shown by figures 3 and 4. 

A sphere represents the best form for which streamlines and 
water-drop trajectories are now known (reference 4) and which would 
have approximately the same rate of water impingements as the front 
of the fuselage. For calctilations a sphere was selected with a 
radius equal to the maximum radius of the fuselage. Uhe relative 
location of the flush test panel on this sphere is shown in. figure 31 . 

Water intercepted by the sphere .- The projected area of the 
sphere is used with the impingement efficiency to determine the weight 
rate of water interception: 

VoTtC^ 

w = 0.225 (19) 

T nn 


where 

C the radius of the sphere, feet 

Weight rate of water intercepted per unit of impingement, area .— 
The weight rate of water intercepted by the sphere divided by the ^ 
surface area over which the water droplets impinge gives the average 
rate of impingement over the area of impingement. 

The equation for the impingement area is : 

Ai = 2rtC2(l - cos ©m) (20 ) 


where ' 

half the central an^le of the total area of impingement on a 
spherical surface 
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Hieref ore, the weight rate of water intercepted per unit of impinge- 
ment area may be written: 


M = 0.225 — Uq m 




100 2rtC^ (1 — cos 0fij) 

which may be reduced to the form. 


M = 0;225 


Tj Uq m 


3-00"^l - cos 9 m) 


( 21 ) 


Weight rate of vater intercepted l)y the flush test panel ,— 

The assmption was made that. Tor watei>-drop sizes of sufficient 
diameter to cause the area of impingement to include the flush test 
panel, the water intercepted per unit area of the panel would he 
approximately equal to the average weight rate of water intercepted 
over the area of impingement. On the basis of this assumption, the 
following equation was written: • 


w = 0.225 


T] Up m A 

i00"^(l - cos 6^) 


( 22 ) 


where 

-A the area of the flush test panel, l.l6 square feet 

Calculated values of water impingement on the flush test panel 
for various water-drop sizes are presented in table YI. 



APPETTOIX B 


Calculations of the Heat Reqviirements on the Flat— Plate 
and Flush Test Windshields for Icing Condition 6 

I. Heat requirement for the flat— -plate vindshield .— 

A, The follcwing assumptions were made: 

1. The water intercepted hy the windshield is the same as that 

intercepted hy a rihhon with an area ec[ual to the 
, projected area of the windshield. (See reference 4 and 
. figure. 30.) 

2. The water that is intercepted by the windshield panel is 

heated to the surface temperatiire and then all, or in 
part, evaporated. ■ . 

3. The windshield surface is ccmpletely wetted. 

4. There is a region of air stagnation at the leading edge 

of the panel. 

Ihe air flow over the windshield is turbulent, 

B. The conditions for illustration are the same as icing condi- 

tion 6, table I: 

, 1. Pressure altitude: 12,700ft 

2. Airplane true airspeed: lUO mph 

3* Free-stream, air temperature; — 2° F 
U. Liquid water content: 1.0 ^am/cu meter 

5. Water-drop diameter: 44 microns 

6. Drop-size distribution: ’ uniform 

, 7. Heated area of panel: 1.49 sq ft 

8, Convective heat— transfer coefficients based on test data 
taken at the center of the panel and the average value 
of s = 1.68 ft. 
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9. Panel angle- vith the fuselage: 45° 

10, Mean surface temperature of the panel: 36° F 

C. External convective heat— transf er coefficient. Since the 

local velocity was measured only at the center of the flat- 
plate panel and eqtiation (15) (Discussion) gives the 
experimental convective heat— transf er coefficient at the 
center of the panel, a conversion factor was needed to 
convert the coefficient at the center of the panel to an 
average value for the entire heated area. By substituting 
2s for I in equation (7) and eqmting the two equations, 
the average coefficient was found to he 1.1 times the 
coefficient at the center of the panel. Therefore, hy using 
equation (15) For this case: 

h = 1.1 X 0.46 (Ui7)°‘® 

T„ = + 460 = 477° F absolute 

*av 2 

Ux = 106 X 1.467 = 156 ft/sec 


7 = 0.0525 Ib/cu ft 

theref ore 


h =17.3 Btu/hr, sq ft, °F 
D. Bate of water impingement, equation (10) : 

T 1 Ap 

M . 0.225 U„ m J 


The percentage of water intercepted is obtained from cxirves 
presented in reference 4 after values of K and q> have been 
calculated from the equations of that reference : 


K 


cp = 


gyya Up 
9iJgb 
l87^Uob 

VS 


= 1.21 


= 1.48 X 10 
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vher e 

7^ specific weight of the water drops, 1 gram/cu cm 

7 specific weight of the air, 8.22 x 10 ^ grams/cu cm 

a radius of the water drop, 0.0022 cm 

b projected height of .windshield or half of ribbon width, 

1.5 X 30.1)3 X cos 1^5° = 32.3 cm 

Uq airplane speed., 6,25P cm/sec 

g acceleration due -to -gravity, 9^0 cm/see, sec 

H absolute viscosity of air, I.69 x 10 gram sec/sq cm 

For these values of K and 9 the impingement efficiency is 
equal to 43 percent, and the equation for the rate of water 
impingement becomes : , 


1 4o X co^ 

M = 0.225 X .43 X 205 X X 1 = 14 Ib/hr, sq ft 

149 

w = MA = 20.9 lb /hr 

E, Calculation of, the factor X, equation (l4) ; 

T = 1 + (ss-Qq) 

For the sample conditions 
Lg = 1074 Btu/lb at tg = 36° F 
Cp = 0.240 Btu/lb, °F 

= 18.55 in. Hg 

e_ = 0.212 in. Hg 

e^ = 0.040 in. Hg 
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t3=-36°F 
= - 2 ° F 
Therefore 

X = 1 + 0.675 = 1.675 

F. Heat required at the outer surface of the flat-opiate panel for 

ice prevention, equation (I3) : 

q = hX(t3-to)+M(t3-to)-0.832 h r(^ 

= 17.3 X 1.675 X 38 + 14 X 38 - 0.832 17.3 X 0.89 (2.05)? 

V 

= 1580 Btu/hr, sq ft , 

G. Temperature of the panel sturface required to evaporate all the 

water intercepted: 

Equating the equation of mass transfer (equation ( 8 ) without 
the latent heat of vaporization) to the weight rate of water 
inipingement, the difference in vapor pressures required to 
evaporate all the water intercepted is first determined: 



As a first approxiiaation, h can he taken as the previously 
calculated value of 17.3 Btu/hr, sq ft, °F . 

then 

eg-e^ = 5.8 in. Hg 

eg = 5.84 in. Hg = 148.3 mm Hg (satirration) 

therefore, 

tg = 139° F 
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For this value of tg, h is recalcvaated to he I7.8 Btu/hr. 

3(1 ft. Op 

therefore, 

eg-©o .= 5.97 in. Hg 

e^ = 6.01 in. Eg = 153 mm Hg 

t = .lkl° -F 

s 

H. Heat required at the outer stirface of the panel to evaporate 
all the water intercepted with no water runhack or blow— off. 

h =17.8 Btu/hr, sq ft, °F 

M =14.0 Ih/hr, sq ft 

tg = 141° F 

X = 7.0 

q = 17.8 X 7.0 X l4l + 14.0 X l4l - 56 = 19,400 Btu/hr, sq ft 

II. Heat requirement for a flush-type windshield . ~ The method 
of calculation and the assumptions made were similar to those used for 
the flat-plate windshield. Icing condition 6, table T, is again 
selected for illustration. For the flush— type windshield the 
assumption was made that the water intercepted by the windshield is 
equal to the water intercepted by a portion of a sphere with a 
radius eqxial to the maximum radius of the upper half of the fuselage . 
(See Appendix A, and fig. 3I.) 


A. External convective heat— transf er coefficient. Since the 

flush windshield test ijanel did not deviate from the fuselage 
contours, the distance from the stagnation region at the 
nose of the fuselage to the center of the panel' was utilized 
in equation' (6), to determine the' convective heat— transfer 
coefficient. 


h = 0.51 T. 


0.3 


av 


s 


Ui7 


0.8 


= 475 F, absolute 


Ui = 190 ft/sec. 

y —0.0525 Ib/cu ft 
= 5.0 ft 


8 
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Therefore, 

h = llv.8 Btu/hr., 8<i ft, °F 

B, Bate of water impingement, equation (22), Appendix A, The 

percentage of water intercepted was obtained from curves 
presented in rsfai’sace II- after the following values of K 
and cp had "been calculated: 

K = 0.2l^2 

cp = 7.5 X 10^ 

then 

T) 5= 3 percent 
and 

The central angle % = l6° does not represent an impinge- 
ment area extensive enough to include all the area of the 
flush windshield (fig. 31 ) > yet test observations indicated 
that during icing condition 6 ice tended to form over the 
entire windshield surface. For this reason, fig = 27.5° 

Wa,s substituted for % in equation (22) . The weight rate 
of water impingement was determined by the equation: 

qUo^ , 

M = ^ = 0.225 r = lb /hr, sq ft 

A jlOO'^1 -cos ©b) 

C. Calculation of -the factor X, equation (ik) 


X = 1 + 


0.622Le 
° p ^o 



+ 0.633 = 1.633 


tg = 33° F 
V - -2°P 
Lg - 1075 Btu/lb 


Pq = 18.55 in. Hg. 



Cp = 0.24 Btu/ll5, °F 
6g = 0,1878 in. Eg 
= 0.04 in. Hg . 

Total heat required at the outer surface, equation (13) . 

- (>-/J 

q = 14.8 X 1.633 X 35 X 35 - O.832 x 14.8 x 0.8(2.05)2 

/o/p 

= ^fe-Btu/hr, sq ft 


LIIRART 
* Of floe of 

^•FOBsatioal lotoiiig^oo 
Batlonal 4dvii^ory Coimltto# 
for Aoronautlco 
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TABLE VI.- CALCULATED VALUES OF WATER IMPINGEMENT 
ON A FIUSH WINDSHIELD PANEL. BASED ON DATA 
CALCULATED FOR SPHERES, REFERENCE 4 


Drop 

diameter 

(cjicrons) 

n 

(percent) 


Sphere 
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(lb /hr) 

Flush 

test 

panel, 

M 

(Ib/h^ 
sq ft) 

Flush 

test 

panel, 

w 

(Ib/hr) 

20 
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0° 
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0 

0 

30 

<1 

3° 

48 



50 

3 

16 ° 




70 
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CO 

OJ 




100 

15 




11^1 


Pressure altitude, 10,000 ft 
True airspeed, I 50 mph . 
Ambient-air temperature, 0° F . 
cp = lil X 105 . 

0M = half of the central angle 
of the total, area of impingemant 
on a spherical surface 
Liquid water content, 1.0 gm/m® ' 


• MTIONAL ADVIS0-.7Y 
COMMITTEE F0''3 A.E70NAUTICS 





























TABLE VII.- COMPARISON OF CALCULATED AND MEASURED . 
HEAT FLOW FROM THE SURFACE OF THE FLAT-PLATE 
PANEL. C-U6 AIRPLANE, WINTER 19k3-h6 
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.TABLE VXII.- COMPARISON OF CALCULATED AND MEASURED HEAT FLOW FROM THE SURFACE OF 
THE V-TYPE WINDSHIELD, C-46 AIRPLANE. WINTER 1946-4-T. 
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TABLE IX,— A CALCT3LATI0 N OF IHE AMOUNT OF WATEE 
intercepted by A FLAT-PLATE PANEL USING TEE 
DROP-SIZE DISTRIBUTION. ICING CONDITION 12 


Liquid 

water 

content 

(percent) 

Liqidd 

water 

content 

(gm/m3) 

Distri- 
bution 
E . 

ax’icf. 
(cm) • 

(percent) 

w 

(Ib/hr) 

5 

0.07 

■ 0,23 

■2.19 

0 

0 

10 

.lU . 

.10+ 

I+.18 

0 

0 

20 

.28' 

.65 

;6.l8 

P , 

0 

30 . 

.h 2 

1.00 

: 9.5'5 

,2 

0,571 

20 

.28| 

i.i+B ' 

14.08 

13 

2.47 

10 


2.00 

19.0 ; 

30 • 

2.855 

5 

.07: ’ 

2.71 

25.8' 

, 

2.24 

100 

1.1+0 

■ : . 

r 

•' r- - 1 

8.136 


4 

cp .= 2.66 X 10 ■ , 

^sed on mean— ef f ectire drop diameter of 19 microns 
ao = 0.00095 cm 
w = 1.91 Ib/hr 
Cl = 890 Btu/hr, sq f t 

Based on drop-size distribution E 
w = 8.l4 lb /hr 
q = 1035 Btu/hr, sq ft 

Experimental data 
q = 1235 Btu/hr, sq ft 

Water :drop^ize distribution and percentages of 

liquid vater content vere obtained from reference h. 
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TABLE X.- CALCULATED MINIMUM WINDSHIELD HEAT REQUIREMENTS FOR A 
MODERATE ICING CONDITION IN TWO TYRES OF ICING CLOUDS;, 
BASED ON AN ASSUMED DROP-SIZE DISTRIBUTION E 
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Figure 1,- The C-46 airplane as equipped for icing research flights in the winter 1945-46, 
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Figure 2.- The C-46 airplane as equipped for icing research flights in the winter 1946-47. 
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Figure 3.- Adjustable flat-plate windshield and flush windshield panel 
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Figure 4.- Electrically heated flush windshield panel, 
viewed from the exterior of the C-46 airplane cockpit. 
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(b) Top view. 

Figure 5.- V-type electrically heated windshield 
installed on the C-46 airplane. 
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oo- Pilot's flush -type' winoshiqlo 

WlMTEe 1^45-46 



FLAT Plate pansi- 
vJiNTee. \aAS-4.® 



CO- PU-OT^ HALF V'TYPe VA> (NO S H\ELO t winter 



FlGUeC 6 - 5KETCH SHOWING THE THER.IAO COUPLE LOCJ¥TIO»^S ON 

VARIOUS u»(N03Hieud configurations use o. on the c-<vg airplane 
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Figure 7.- Single panel from the V-type windshield installation. Heating 
was provided by transparent^ electrical— conducting film under outer 
glass layer. 
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Figure 8.- Details of external discharge heated-air ice-prevention system for 
flush windshield of C-46 airplane. 
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TBUE: ftie >SPE£0 
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fa) In icing conditions. Time, 1549. 



(b) Immediately after leaving icing conditions. Time, 1551. 


Figure 13.- Flush windshield of C-46 airplane during icing 
condition 6, table I, with inadequate heat supply to the 
external discharge system for the pilot ^s panel. 
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TIME 

PRESSURE ALTITUDE 

TRUE AIR SPEED 

AMBIENT AIR TEMPERATURE — 

LIQUID WATER CONTENT 

MEAN EFFECTIVE DROP SIZE- 

HEAT SUPPLIED TO PANEL 

HEAT SUPPLIED TO WINDSHIELD 


1650 TO 1720 

10,300 TO 10,600 FEET 

170 TO 150 MPH 

20 TO I7.5*F 

.15 TO .3 GRAM/M* 

50 TO 20 MICRONS 

1205 BTU/HR FT^ 

525 BTU/HR FT* 


Figure 14.- Sketch based on photographs and flight engineer’s notes 
of ice accretions on nose and windshields of C-46 airplane after 
icing condition 9, table I. Panel angle 60° with the fuselage. 


Page intentionally left blank 


Page intentionally left blank 



NACA TN No. 1434 


77 



(a) In icing conditions, start of removal. Time, 1703. 



(b) In icing conditions, removal continued. Time, 1706. 

Figure 15.- Windshield of C-46 airplane during icing condition 9, 
table I, showing removal of ice accretions from the electrically 
heated flush panel, right side, and the external discharge area, 
left side. 
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(c) In clear air, removal completed. Time, 1720. 



(d) Ice formation on ^-inch-diameter rod during the 
icing run. 


Figure 15.- Concluded. 
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Figure 16.- Windshield of C-46 airplane after a flight in 
icing conditions in which there was no tendency for ice 
to form on the windshield. Ice accretions shown formed 

on --inch-diameter rods beside windshields; icing 
8 

conditions 10 and 11, table I# 
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FLIGHT NUMBER 31 
PRESSURE ALTITUDE 
TRUE AIR SPEED 
AMBIENT AIR TEMPERATURE 
LIQUID WATER CONTENT 
MEAN DROP SIZE 
HEAT SUPPLIED TO PANEL 
ICING CONDITION NUMBER 


TIME 1116 
8,000 FEET 
160 MPH 
14 *F 
.2 GRAM/M3 
10 MICRONS 
1520 BTU/HR FT^ 




FLIGHT NUMBER 31 
PRESSURE ALTITUDE 
TRUE AIR SPEED 
AMBIENT AIR TEMPERATURE 
LIQUID WATER CONTENT 
MEAN DROP SIZE 
HEAT SUPPLIED TO PANEL 
ICING CONDITION NUMBER 


TIME 1628 
11,650 FEET 
180 MPH 
17 -F 
.5 GRAM/M3 
18 MICRONS 
1240 BTU/HR FT* 

2 




FLIGHT NUMBER 33 
PRESSURE ALTITUDE 
TRUE AIR SPEED 
AMBIENT AIR TEMPERATURE 
LIQUID WATER CONTENT 
MEAN DROP SIZE 
HEAT SUPPLIED TO PANEL 
ICING CONDITION NUMBER 


TIME 1147 
8,300 FEET 
170 MPH 
18 *F 
.9 GRAM/M3 
10 MICRONS 
1230 BTU/HR FT2 
3 NACA 

A-i 0 8 I 


Figure 17.^ Ice accretions formed on the flat-plate windshield of the 
C-46 airplane during icing conditions 1, 2, and 3, table !• Panel 
angle 30^ with the fuselage. 
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FLIGHT NUMBER 2 3 
PRESSURE ALTITUDE 
TRUE AIR SPEED 
AMBIENT AIR TEMPERATURE 
LIQUID WATER CONTENT 
MAXIMUM DROP SIZE 
HEAT SUPPLIED TO PANEL 
ICING CONDITION NUMBER 


TIME 1454 
10,350 FEET 
135 MPH 

.6 TO .B GRAM/M3 
II TO— MICRONS 
910 BTU/HR FT* 




FLIGHT NUMBER 23 
PRESSURE ALTITUDE 
TRUE AIR SPEED 
AMBIENT AIR TEMPERATURE 
LIQUID WATER CONTENT 
MAXIMUM DROP SIZE 
HEAT SUPPLIED TO PANEL 
ICING CONDITION NUMBER 


TIME 1516 
11,500 FEET 
145 MPH 
-0.5 *F 
.3 TO 1.2 GRAM/M3 
20 TO — MICRONS 
16 30 BTU/HR FT* 
5 




FLIGHT NUMBER 23 TIME 1548 

PRESSURE ALTITUDE 12,700 FEET 

TRUE AIR SPEED 140 MPH 

AMBIENT AIR TEMPERATURE -2 »F 

LIQUID WATER CONTENT .6 TO 1.4 6RAM/M3 

MEAN DROP SIZE IS TO 44 MICRONS 

HEAT SUPPLIED TO PANEL 1720 BTU/HR FT* 
ICING CONDITION NUMBER 6 naca 

A 10 982 

»2-» 1-46 


Figure 18.- Ice accretions formed on the flat-plate windshield of the 
C-46 airplane during icing conditions 4, 5, and 6, table I. Panel 
angle 45^ with the fuselage. 
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FLIGHT NUMBER 49 
PRESSURE ALTITUDE 
TRUE AIR SPEED 
AMBIENT AIR TEMPERATURE 
LIQUID WATER CONTENT 
MEAN DROP SIZE 
HEAT SUPPLIED TO PANEL 
ICING CONDITION NUMBER 10 


TIME 1121 
9,140 FEET 
165 MPH 
9 5 *F 
.7 GRAM/M3 
23 MICRONS 
1170 BTU/HR FT* 




FLIGHT NUMBER 49 TIME 1132 

PRESSURE ALTITUDE 8,900 FEET 

TRUE AIR SPEED 165 MPH 

AMBIENT AIR TEMPERATURE II 

LIQUID WATER CONTENT .2 GRAM/M3 

MEAN DROP SIZE 17 TO 44 MICRONS 

HEAT SUPPLIED TO PANEL 1190 BTU/HR FT* 
ICING CONDITION NUMBER II 




FLIGHT NUMBER 49 
PRESSURE ALTITUDE 
TRUE AIR SPEED 
AMBIENT AIR TEMPERATURE 
LIQUID WATER CONTENT 
MEAN DROP SIZE 
HEAT SUPPLIED TO PANEL 
ICING CONDITION NUMBER 12 


TIME 1155 
10,600 FEET 
160 MPH 
5 •F 
1.4 GRAM/M3 
19 MICRONS 
23 5 BTU/HR FT* 

NACA 
A 10 9 8 3 
12-1 1-46 


Figure 19.- Ice accretions formed on the flat-plate windshield of the 
C-46 airplane during icing conditions 10, 11, and 12, table I. 

Panel angle 60® with the fuselage. 
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Figure 20.- V-type windshield of C-46 aii 7 >lane during flight 
in icing conditions showing unsuccessful ice-removal 
operation of external discharge jets over pilot's panel, 
with heat supply of 18,000 Btu per hour. 
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Figure 21.- Typical ice accretion on V-type windshield 
without any heat supplied. 



Figure 22.- Ice accretions formed on V-type windshield of 
the C-46 airplane during icing conditions 15 and 16. 
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(a) In icing conditions. Time, 1140. 



(b) In clear air, ice removed. Time, 1215. 

Figure 23.- Ice accretions fomed on V-type windshield of the 
C-46 airplane during icing condition 17. 
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FIGURE 25 -TEIvnPERaTU BE GRADIENTS IN THE HEATED BOUNDARY 
LAYER OVERTHE FLUSH TYPE WINDSHIELD OF THE C-46 A-IRPlANE. 
PEESSURE ALTITUDE )2>000 FEET ,TRUE AlR SPEED IS'S’ KPH. 
HEATED AlP SUPP U ED 13,000 BTU/ HE WINTER »a45'-4fe 
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